LEAF LITTER DECOMPOSITION ESTIMATES OF GLOBAL 
VARIABILITY BASED ON YASSO07 MODEL 



M. TUOMIi'2, T. THUM3, H. JARVINEN^, S. FRONZERi, B. BERG-i'5, M. HARMON^, 
J. A. TROFYMOW'', S. SEVANTO*, AND J. LISKI^ 



Abstract. Litter decomposition is an important process in the global carbon 
cycle. It accounts for most of the heterotrophic soil respiration and results in 
formation of more stable soil organic carbon (SOC) which is the largest ter- 
restrial carbon stock. Litter decomposition may induce remarkable feedbacks 
to climate change because it is a climate-dependent process. To investigate 
the global patterns of litter decomposition, we developed a description of this 
process and tested the validity of this description using a large set of foliar 
litter mass loss measurements (nearly 10000 data points derived from approx- 
imately 70000 litter bags). We applied the Markov chain Monte Carlo method 
to estimate uncertainty in the parameter values and results of our model called 
Yasso07. The model appeared globally applicable. It estimated the effects of 
litter type (plant species) and climate on mass loss with little systematic er- 
ror over the first 10 decomposition years, using only initial litter chemistry, 
air temperature and precipitation as input variables. Illustrative of the global 
variability in litter mass loss rates, our example calculations showed that a 
typical conifer litter had 68% of its initial mass still remaining after two de- 
composition years in tundra while a deciduous litter had only 15% remaining 
in the tropics. Uncertainty in these estimates, a direct result of the uncertainty 
of the parameter values of the model, varied according to the distribution of 
the litter bag data among climate conditions and ranged from 2% in tundra to 
4% in the tropics. This reliability was adequate to use the model and distin- 
guish the effects of even small differences in litter quality or climate conditions 
on litter decomposition as statistically significant. 
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1. Introduction 



Litter decomposition plays a crucial role in the global carbon cycle. Carbon 
dioxide emissions from the decomposition of soil organic carbon (SOC) are equal 
to about 60 Pg of carbon per year, whi ch is about s even times as much as the 
annual emissions of fossil carbon (IPCC reportl 120071 ). Most of these emissions 
originate from the decomposition of the relatively labile litter. Litter decomposition 
results also in formation of more stable org anic compound s . Th ese compounds 
represent the majority of all SOC (Davidson and Janssens . 2006), which is the 
large st terrestrial carbon stock, equal to about 2300 Pg (IJobbagv and Jackson . 
2000), or three times the atmospheric carbon stock today ( IPCC report . 2007 ). 



Litter decomposition will respond to changes in climate and, because of its great 
importance to the global carbon cycle, there may be remarkable feedbacks to the 
future climate change. Litter decomposition is also considered a complex process 
controlled by numerous other factors. To understand this process and to improve 
the estimates of the role of litter decomposition in the global carbon cycle, a global 
model of litter decomposition is needed. 

A variety of approaches has been applied to estimate litter decomposition at large 
geographical scales. These approaches can be divided into three groups, 1) regres- 
sion models ba sed on lit ter bag s tudies (e.g. Meentemever . 1978t Berg et al. ■ Il993t 
Trofvmow et al... 2 002: Zhang et al. . 2007.) or on soil respiration measurements (e.g. 
Llovd and TavloI^ll994^ ■ 2) sp ecific dynamic soil carb on models comprising of com- 



partments (e.g. iParton et a. 



i pecinc dynamic sou carb on models comprising oi com- 
I llMZl; [J enkinson, 'l99d) or models based on a theory 



of continuous SOC quality ( Bosatta and Agren . 2003, ) and 3) less specific dynamic 
soil carbon models used for investigating the dynamics of nitrogen mineralization 
( Parton et al. I. l2007 l: Manzoni et al., 1 20081) or applied to national carbon account- 



ing (e.g. Kurz and Adds . 2006: Lis ki et al.L [ioOS. 'l or to dynamic glob al vegetation 



mode ls (e.g. ISitch et al.l . l2003l ) or Earth System Models (ESM) (e.g. Ijones et al 
200l 



A particular challenge in developing SOC decomposition models is that the in- 
ternal SOC pools of the models and especially car bon fluxes between the pools can- 
not easily be determined from measurement data (IChris tensen . 1996t [Elliott et all . 
[1996). Therefore, it is a common practice to make a decision of the fluxes flrst and 
then quantify the para meters determining the magnitudes of these fluxes based on 



measurement data (e.g. Moor head et al. . 19991 ). 



To investigate the global patterns of litter decomposition, we wanted to avoid 
this uncertainty stemming from the preflxed model fluxes and assumed that any 
fluxes between the SOC pools of our model were possible. Then, we determined 
these fluxes and the values of all other model parameters directly from a global 
data set of l itter m as s loss measurements ( Gholz et al. [ l2000l: l iTrofvmow et al. [. 119981: 
Berg et al.[ . [l991a[ [bl. [1993*) using the Markov chain Monte Carlo (MCM C) method 



with Metropolis-Hastings algorithm ( Metropolis et al. . 1953 : Hastingsl , [l970 ). We 
also required that the same parameter values should flt to the entire global data set 
so that our litter decomposition model would be applicable across climate conditions 
worldwide and to a wide variety of litter types. 

Our aim was to calculate the full posterior probability density of model parame- 
ters, i.e. the joint probability density of all the free parameters in the model given 
the measurements. Possibilities to calculate such statistical uncertainty estimates 
have been lacking from previous SOC models although they are necessary if we want 
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to evaluate the reliability of model results in a statistical sense. The uncertainty 
estimates can only be calculated in a meaningful way if the model is not over- 
parameterized. When building the model we required that there were no modelling 
errors introduced by over-parameterization. We also required that the model we 
developed fulfilled the Occamian principle of parsimony, which means that we chose 
the simplest model structure from a set of almost equally good alternatives. These 
requirements were made to find a reliable model whose prediction uncertainties are 
known reliably. 

The objectives of this study were, first, to develop a global model of foliage litter 
decomposition that met the above requirements, and second, to estimate the global 
patterns of foliage litter decomposition and their uncertainty using this model. 



2. Material and methods 



2.1. Measurements. To develop the global model of foliage litter decomposition, 
we used litter bag data from 97 sites across Europe, and North and Central America 
(Table [2]). The d ata from the USA and Central America were collected within the 
LIDET n etwork (ICholz et al.l. l2000l) and the data from Canada within the CIDET 
network (jTrofvmow et al.L Il998h. wh e reas th e data from Europe originated from 



several research projects ( Berg et al. . 1991al lb. 1993 ). The study sites covered a 



wide range of climate conditions in terms of temperature and precipitation (Fig. 
[TK), which are the most important c l imate factors affecting l itter decomposition 
(lMeentemevelll978tlB"erg et al.l . [l99llAert£l[l997tlLiski et"al] . l2003l : IParton et all . 

2nn7h . 



The litter bag datasets consisted of measurements for foliage litter of 34 plant 
species including several coniferous and deciduous trees (Table[3]). The initial chem- 
ical composition was meas ured for each litter type and the loss of mass was fol 
lowed fo r 3.1 to 10.2 years (jCholz et al.l . I2OOOI : iTrofvmow et al.l . Il998l: iBerg et al ' 



1991al lbl. [l99l . In addition to the total mass loss, the mass loss of chemical co m- 
pound groups was measured at seven Swedish study sites ( Berg et al. . 1991al |bh. 
Together the datasets comprised 9605 data points. The LIDET data consisted of 
values for individual litter bags, whereas the Canadian data p oints were averages 
of four litter bags collected at the same time ( Trofvmow et al.l 12002). M ost of the 
European data points were averages of 25 litter bags (|Berg et al.l . ll991al lbh. Thus, 
the data was received using approximately 70 000 litter bags. 

In addition to the litter bag data, we used a dataset on accumulatio n of SOC at 



1998 



26 sites alo ng a 5500 year long soil chronosequence in southern Finland (jLiski et al 
2005^). This data provided us with both information on formation of humus 



from decomposing litter and humus decomposition. 

2.2. The model. The Yasso07 mod el developed in th is study is a generalization 



2003). Yasso07 is based on three 



of an earlier Yasso soil carbon model (jLiski et al 
assumptions of litter decomposition: 

(1) Non-woody litter consists of four compound groups, i.e. compounds soluble 
in a non-polar solvent, ethanol or dichloromethane (denoted using E), or in 
water (W), and compounds hydrolysable in acid (A) and neither soluble nor 
hydrolyzable at all (N) . Each group has its own mass loss rate independent 
of the origin of the litter. These compound groups are called the labile 
groups. 
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Figure 1 . a) Distribution of the global cli mate conditions on land 



(black dots) from the CRU 2.1 database (|New et all . 120021 ) and 
the 97 litterbag measurement sites (LIDET sites red, CIDET sites 
blue and EURODECO sites green). Variables T,n and Pa are the 
mean annual temperature and annual precipitation, respectively. 
LIDET, CIDET and EURODECO measurement sites are in the 
USA and Central America, Canada, and Europe, respectively, b) 
95% confidence intervals as a function of climate for estimates of 
mass remaining after two years of decomposition (a unitary initial 
mass) calculated using the YassoO? model. The error is shown for 
100 randomly selected climatic conditions. 



(2) The mass loss rates of the compound groups depend on the climatic condi- 
100 tions that can be described simply by using temperature and precipitation. 

(3) Decomposition of the compound groups results in mass loss from the system 
and in mass flows between the compound groups. In addition, the mass 
loss of the four compound groups results in formation of more recalcitrant 
humus (H). 
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The first assumption is justified by ea rUer studies showi ng; that the above chem- 
ical groups decompose at different rates ( Berg et al] . ll982l ). and that this g^ r oupiii g 
differentiates fitter types according to decomposition rate ( Palosuo et a l.. '20051}^ 



The second assumption was based on r esults of several earfier studies (e .g. Meentemever 



1978: Berg et al 



i )tion was ba; 
1993t lAertd . 



1997t iLiski et all 12003 : iParton et all . 120071 ). The 



third assumption follows from a general view on carbon cycling in soil that in- 
volves transformations of organic compounds. This includes the break-down of 
complex compounds and the formation of simpler ones in chemical decomposition 
reactions, which are catalyzed by enzymes excreted by soil microbes; the uptake of 
the simplest compounds by soil microbes; and the formation of carbon dioxide and 
synthesis of biomass in the metabolism of soil microbes. 

Yasso07 was formulated according to these hypotheses. All parameter values 
were treated as free parameters when the model was fitted to the data. Hence, 
Yasso07 is a set of first order differential equations and defined as 

(2.1) x(i) = A(C)x(t) + h{t), x(G) ^ Xo 

where x = {xa,xw,xe,xn,xh)'^ is a vector describing the masses of the five 
compartments as a function of time (i); A(C) is a matrix describing the de- 
composition rates and the mass fiows between the compartments as a function 
of climatic conditions (C); vector h{t) is the litter input to the soil. Vector 
Xo = ixA,o,xw,o,XE,o^XN,o,XH,o) is the initial state of the system and Xi,o the 
initial chemical composition, with ireferring to A, W, E, N and H. Matrix A is 
defined as a product of the mass flow matrix Ap and the diagonal decomposition 
coefficient matrix k(C) = diag(fcyi, /c^y , fc^;, /cjv, fc// )(C), where ki are the decom- 
position rate coefficients of the compartments. The mass flow matrix is defined 
as 



(2.2) 
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where pi € [0, 1] are the relative mass fiow parameters between the compartments. 
It is further assumed that ki — fci(C), for all i, and that the sum of parameters pi, 
describing the mass flows out of any of the compartments, does not exceed unity. 
This mass fiow matrix was selected because all the possible flows between A, W, 
E and N are present. This is the most general model structure that can still be 
presented as a linear differential equation (Eq. 12. ip . 

The climate dependence of the decomposition rate factors ki is formulated as 



(2.3) 



h{C) = a, exp (/3iT + (i^T^) (l - exp[7Pj), 



where T is temperature (Celcius scale) and Pa is the annual precipitation and a^, 
(3i , (32 and 7 are free paramete r s. Th is form of the temperature dependence was 
justified earlier bv lTuomi erall (|2008l ). 

We also tried to determine the separate values of these free parameters for each 
compound group, but it resulted in an over-parameterized model and a lower pos- 
terior probability. This happened because mass loss was measured by compound 
group only at the Swedish study sites (see Table [3]). Conseqeuently, we had to 
apply the same parameter values for each compound group. 



6 



M. TUOMI ET AL. 



The intra-annual variations in temperature resulting from seasonal changes were 
approximated by using a sinusoid 

(2.4) T(t) = T™ + T„ sin (^), 

130 where Ta = ^{T^ ,max 2^m,rTiin) <ind Tj^^^^nax &nd T^^^^n 9.re the maximum and 
minimum mean monthly temperatures, respectively, and tp]s the period of one 
year. 

The full inverse solution, i.e. the probability density of the parameter vector 6*, 
was found by sampling the parameter space with the MCMC method. After a burn- 
135 in period, the proposal density was constructed to allow for the convergence of the 
chain to the posterior density with an acceptance rate of approximately 0.3. With 
the full probability density available, we then calculated the maximum a posteriori 
(MAP) point estimate and Bayesian 95% confidence set for the model parameter 
vector. 

When calculating the posterior density, we used a Gaussian likelihood function 
for the measurements, such that 

(2.5) f{e\m) (X exp 1^ - i ^ ^,5, j ^{6), 

140 where f{9\m) is the conditional posterior density of the model parameters with 
the condition that measurements m = (rni, ...,mj, ...) have been made and 'k{6) is 
the prior density of the parameters. Variable Si is the common sum of squared 
residuals of measurement set and is some relative weight of measurements 
mi. These weights were set to balance between data sets with large number of 

145 measurements containing little information and data sets with few measurements 
but high information content. 

Finding the MAP estimate is equal to finding that approximately satisfies the 
conditions 

(2.6) 6 = arg min uJiSj 

i 

and 

(2.7) ^^argminS'i 



for all the data sets Wj simultaneously. The weights uji were selected so that 
the latter condition is satisfied as well. This is basically made to extract all the 
information on model parameters from the measurements and to make siirc that 
150 there is as little systematic error between the model and any of the measurement 
sets as possible. 

We chose not to calculate the standard goodness-of-fit measures, the -values 
for the different datasets, because they provide no extra information in addition to 
the residuals between model with MAP parameter values and measurements. Also, 
155 to use the values, they should be compared with corresponding values calculated 
using different models. As there are no other models that predict the climate 
dependent time-evolution of the A, W, E, N and H compounds, this comparison 
could not be made. 

We tested several alternative structures within the model framework set by our 
three assumptions, such as different precipitation dependence functions or different 
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temperature or precipitation dependences for the different compound sets. Tfiese 
structures were compared using the model probabilities according to the Bayesian 
model selection theory, which automatically contains the Occam's razor and there- 
fore penalises unnecessarily complicated model structures. These probabilities were 
defined as 

(2.8) Pig^\n^) = ^MEjEM_ p(^^,\g^)= f f{m\e,, gM0,\9,)d9,, 

where gj is the jth model structure, 9j G f2j are its parameters, P{gj\m) is its 
likelihood and P{gj) its prior probability. These prior probabilities were all set 
equal. We simply selected the structure that has the highest probability accord- 
ing to Eq. p.Sp with respect to all the measurements available. Desp ite that 



Bayesian methods have proven usefull i n sev eral statistical problems (e.g. lEllison 



2004 iLichter et all . l2005t iTuomi et all . |2008|) . according to our knowledge, this is 



the first time Bayesian methods are applied to the decomposition process of foliage 
litter. 

When fitting the parameters of YassoO? model to the data, we noticed that 
the level of mass loss rate was higher in the European litter bag measurements 
compared to the North o r Central American ones. The same difference has been 



observed already earlier (jPalosuo et al.l . l2005f ) . It is probably caused by a larger 



mesh size and smaller litter mass used in the European litter bags. This may have 
caused more leaching and thus overestimates of mass loss in Europe. To account for 
this difference, we introduced a scaling factor to the model and determined its value 
together with the other parameter values. This scaling factor had a narrow normal 
distribution with a mean equal to 0.58 and 95% of the probability distribution 
between 0.56 and 0.60. The need to include the scaling factor in the analysis adds 
uncertainty to the results of this study. However, this uncertainty is related mostly 
to the estimated level of decomposition, whereas the estimated effects of climate 
and litter quality on decomposition are less dependent on the litter bag type. This 
happens because there were no systematic differences in the climate or litter quality 
effects between Europe and North and Central America despite the differences in 
the litter bags. 

3. Results 

We were able to establish a global description of foliage litter decomposition pro- 
cess using the large data set of litter bag measurements and advanced mathematical 
methods of complex inverse problems. 

The model fitted the data with little systematic error with respect to any variable 
investigated. The residuals, i.e. the differences between model-calculated estimates 
and measurements, did not deviate significantly from zero for 31 litter types out of 
34 when looking at the 68% confidence intervals (Fig. [2^). With a 95% confidence 
level, none of the litter types deviated from zero (not shown) For one of the deviating 
types {Kobresia myosuroides), there were only 20 data points, which makes it more 
probable that these measurements deviated by pure chance. The residuals were 
not correlated with either the time elapsed since the start of decomposition or 
any climate variable (Fig. [2lD-e). The residuals were also uncorrelated with the 
initial nitrogen concentration of the litter species despite the fact that we did not 
include the effect of nitrogen in the model (Fig. [2f). The nitrogen concentration 
ranged from 0.3 to 2.2%. The lack of systematic error means that the description of 
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Figure 2. Model residuals (Data - Model) by litter type (a), linear 
1st order systematic trends in the residuals by litter type as a 
function of time since the start of decomposition (b), annual mean 
temperature (c), annual precipitation (d), difference between the 
coldest and warmest months of the year (e) , and the residuals as 
a function of the initial nitrogen concentration of the litter types 
(f). The error bars are la values, and thus they cover 68% of the 
probability density. Letters E, C and L after dataset names refer 
to EURODECO, CIDET and LIDET, respectively. 



foliage litter decomposition in the model is, on average, valid for the global climate 
200 conditions and the range of litter characteristics covered by the litter bag data. The 
error in the results of the model is thus of random nature which can be characterized 
by probability densities with zero mean. 
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Table 1 . Maximum a posteriori parameter values of YassoO? and 
the 95% confidence limits. 



Parameter 


Value 




Unit 


Interpretation 




0.66±0.11 






-1 


decomposition rate parameter of A 


^ vv 








-1 


decomposition rate parameter of W 


OiE 


0.35±0.08 






-1 


decomposition rate parameter of E 


OiN 


0.22±0.06 




a" 


-1 


decomposition rate parameter of N 


Pi 


0.32±0.08 








relative mass flow magnitude, W ^ A 


P2 


o.oi+°;t 








relative mass flow magnitude, E — > A 


P3 


0.93+°?? 

— U.ii 








relative mass flow magnitude, N — s- A 


TiA 


34+° l? 








relative mass flow magnitude, A — > W 


P5 


0.00+°°^ 

— u.Uu 








relative mass flow magnitude, E W 


P6 


00+° °^ 

'^•^"-0.00 








relative mass flow magnitude, N W 


P7 


00+° °^ 








relative mass flow magnitude, A E 


PS 


'^■^"-0.00 








relative mass flow magnitude, W E 


P9 


01+° "^ 

'^■"-'^-0.01 








relative mass flow magnitude, N ^ E 


PlO 


00+" °^ 








relative mass flow magnitude, A ^ N 


Pn 










relative mass flow magnitude, W ^ N 


Pl2 










relative mass flow magnitude, E — > N 


Px 


7.6±2.0 


10- 


-2 




temperature dependence parameter 




-8.9±6.5 


10- 


-4 


0^-2 


temperature dependence parameter 


7 


-1.27±0.20 




m 


-1 


precipitation dependence parameter 


PH 


0.040±0.009 








mass flow to humus 


an 


q q+0.6 
0.0_Q 7 


10 


1-3 


a-1 


humus decomposition coefficient 



Fitting the model to the data revealed that there were four major carbon fluxes 
between the labile compartments (Table [T]). Two of them were towards the com- 

205 partments with higher decomposition rates representing breakdown of complex 
compounds and formation of simpler and labile ones during the decomposition 
process. The other two were towards the other direction indicating a re-synthesis 
of more complex compounds in the process. The other fluxes between the labile 
compartments were negligible. Four per cent of the combined mass loss of the labile 

210 compartments resulted in humus formation. 

The decomposition rates were significantly different between the five compound 
groups included in the model (Table [T|). The decomposition rate of the water 
soluble compounds was an order of magnitude higher than the corresponding rates 
of the other three labile compound groups. Among these groups, the decomposition 

215 rates decreased from the acid hydrolyzable compounds to the non-polar extractable 
and the acid unhydrolyzable compounds. The decomposition rate of the humus 
compartment was two orders of magnitude lower. 

The calculated estimates for the mass-loss rates of foliage litter were statistically 
significantly different among sites selected to represent tundra, the boreal zone, the 

220 temperate zone and the tropics and, within each site, among a typical coniferous 
and deciduous litter types (Fig. [3]). The probability densities of these estimates 
represent uncertainty caused by uncertainty about the parameter values of the 
model. The uncertainties of model predictions in Fig. [T|d show the magnitude of 
95% confidence intervals as a function of climatic conditions for a decomposing 
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Figure 3. Probability densities of YassoO? estimates for litter 
mass-remaining after two years of decomposition (a unitary ini- 
tial mass) in different climate zones. A characteristic coniferous 
and deciduous species selected for each climate zone, tundra Picea 
mariana (PM) and Betula papyrifera (BP), the boreal zone Pinus 
banksiana (PB) and Populus tremuloides (PT), the temperate zone 
Thuja plicata (TP) and Fagus grandifolia (FG), and the tropics 
Pinus elliottii (PE) and Drypetes glauca (DG). The latitude, lon- 
gitude, T,n, Ta and Pa of these sites are: tundra 56.32°, -94.85°, 
-4.0°C, 10.1°C, 0.50m; the boreal zone 53.22°, -105.97°, -2.4°C, 
9.7°C, 0.59m; the temperate zone 35.00°, -83.50° 13.9°C 5.1°C, 
1.72m; the tropics 9.17°, -79.85°, 26.0°C, 0.6°C, 2.86m. 

225 body with an unitary initial mass. This uncertainty followed the distribution of the 
litter bag data across the global climate conditions (Fig. [T]) . It was smallest for the 
boreal and temperate zones where most of the study sites were located, and largest 
for the coldest tundra and the tropics which had only a few study sites. 

Illustrative of the global variability in the mass-loss rate of foliage litter, the 

230 coniferous litter of tundra {Picea mariana needles) had about 68% of its mass still 
remaining after two years of decomposition, while the deciduous litter of the tropics 
{Drypetes glauca leaves) had only about 15% remaining (Fig. [3]). The values of 
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Figure 4. Estimated geographic variations in the leaf htter de- 
composition rate. A unitary rate corresponds to chmatic condi- 
tions where — 0°C, Ta = 0°C and Pa is sufficiently large to not 
limit the decomposition. 

the decomposition rate parameters of the model were up to 13 times higher in the 
tropics compared to the northern boreal zone (Fig. [4]). The decomposition rates 
235 arc not directly comparable to the mass loss rates because decomposition results 
also in transfer of carbon between the compartments inside the model. 

4. Discussion and conclusions 

The aim of this study was to estimate the global patterns of foliage litter decom- 
position after developing an appropriate model for this purpose. We required that 
240 the model was globally applicable and the results were accompanied by uncertainty 
estimates. 

We found that the results of the YassoO? model did not deviate systematically 
from the mass loss measurements of different litter types across the global climate 
conditions (Fig. ^ for 31 litter types out of 34 when looking at the 68% confidence 

245 intervals (Ict) and for none of the litter types when looking at the 95% intervals 
{2a). The lack of the systematic deviation means that the probability densities of 
the results of the model represent the random uncertainty caused by uncertainty 
about the parameter values of the model (Fig. [T]3, Based on these results, we 
conclude that YassoO? met the requirements we set for the global model. 

250 We formulated the YassoO? model according to three assumptions of foliage 
litter decomposition (see The Model). Comparison of the model results to the 
measurements provides us with means to evaluate the validity of these assumptions. 

Regarding the first assumption, the decomposition rates of the four labile com- 
pound groups differed statistically significantly from each other (Table [1]). When 

255 applying the same group-specific rates independent of litter type, the YassoO? model 
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estimated the mass loss rate of 31 litter types out of 34 without a statistically sig- 
nificant systematic error (Fig. [2K)- These results support the validity of the first 
assumption. In addition, they suggest that the mass loss of different litter types 
can be estimated based on the initial division of the mass to the four compound 
groups over, at least, the first 10 years of decomposition, which was the time period 
covered by our data. 

Regarding the second assumption, temperature and precipitation appeared in- 
deed as the major climate factors affecting foliage litter decomposition at the global 
scale. After their effects were accounted for, the results of the YassoOT model did 
not deviate from the mass loss measurements in any climate-related way (Fig. [2)3 
to e). Temperature and precipitation, combined in various ways, ha ve been used to 
expla i n the effect s of climate on litter dec omposit i on alr e ady ear lier ( McentcmcyeJ, 
19781 : iBerg et all . Il993l : lAertd . [1997: Liski et all . l2003t IParton e t al., 1007). We 



tested also different combinations of these two variables to find the final and best 
model. The Gaussian formulation of the temperature effects was based on our ear- 
lier comparison of alternative formulations ( Tuomi et al.. . ,2008i) . Unfortunately, the 
data did not allow us to test whether the effects of temperature or precipitation dif- 
fered between the compartments. Trying to fit separate values of these parameters 
resulted in an over-parameterized model and a lower posterior probability. 

Regarding the third assumption, the information in the measurement, based on 
the parameter probability densities, was sufficient for distinguishing four major car- 
bon fluxes between the labile compartments and a fiux to the humus compartment 
in addition to the fluxes out of the system. These fluxes were identified without 
direct measurements of the fiuxes themselves. Yet, our n iodel fits to a general view 
of microbial litter decomposition and humus formation (jStevenso n|, il982 ). Two of 
the internal fiuxes were towards more labile compartments representing probably 
a break-down of complex compounds by exoenzymes exudated by soil microbes. 
The fiux from the water soluble materials to acid hydrolysable ones may represent 
synthesis of more complex compounds in the metabolism of the soil microbes. The 
flux from the ethanol soluble compartment to the non-soluble one is more difficult 
to explain. It may be a result of the fitting procedure of the model to the data 
that led to rather similar decomposition rates of these compartments. Recalcitra nt 



humus is formed from decomposition products in chemical reactions (|Stevensonl . 
[1982) explaining probably the flux to the humus compartment. 

Process-modellers may consider it as a weakness that we applied the strongly 
empirical method without a priori information to determine the fluxes of the model. 
However, it was the only method that we could use because direct mea surements of 
the flu xes are lacking, especially at this extensive geographical scale. Ed air et al 



([20081) ignored fluxes between carbon pools in their model when studying litter 
decomposition across North a nd Centra l Ame rica. However, the common view 
of the decomposition process (|Stevensonl . Il982h and more detailed decomposition 



models indicate that such internal fluxes exist (jParton et al.l . 119871 ) 



Determining fluxes between compartments is a common problem in developing 
dynamic comp artmental soil carbon models because these fluxes cannot usually be 
measured (e.g. ChristensenL 1996 : Elliott et al. . 19961 ). For this reason, it is common 
to d ecide about the fluxes before determining the parameter values of the model 
(e.g. Moorhead et al. ,[i99i). This practise adds uncertainty to the structure of the 
model that is difficult to control. To avoid this problem, we determined each fiux 
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of the YassoOT model from the measurements (Table [3]) . The ability to indentify 
the fluxes from the rather unspecific litter bag measurements is an indication of the 
power of the mathematical methods we used. 

It was surprising that the YassoO? model fitted to the data as well as it did 
even though the nitrogen concentration of the litter types or other characteristics 
of nitrogen availability at the sites, such as site fertility, were not accounted for 
in the model at all. Moreover, the residuals were not correlated with the initial 
nitrogen concentration of the litter species (Fig. [2f). Nitrogen is considered as 
an important regulator of deco mposition and in many other models it is used as a 
controlling input variable (e.g. iParton et all . llQSTt [jenkinson . 19901 ). According to 
the results of this study, climate and chemical litter quality are still stronger controls 
of decomposition at the global scale. Once their effects were modelled carefully, the 
effects of nitrogen could not be distinguished from the remaining variability in the 
data. Nitrogen effects may still be important at smaller scales or in more detailed 
experiments. 

It is difficult to compare the reliability of litter decomposition estimates calcu- 
lated using YassoOT and other models. The reason is that error in the results of the 
other models has not been assessed in a similar probabilistic sense. The regression 
models of litter decomposition have bee n evaluated mainly on the ba sis of Pear- 

an d littl e 
2008h . 



son correlation values or values (e.g. iBerg et al.l . Il993t lAerta. 119971) 



attention has been paid to the reliability of the results (e.g. iTuomi et al 



325 Even if some values were available for other globally applicable models fitted 
to similar measurements, they could not be used for comparison easily. This is 
because the values depend on error estimates of individual measurements and if 
these errors have not been measured they need to be approximated or estimated by 
the modeller. More importantly, however, the values tell us nothing about the 

330 possible over-parameterization of the model. Increasing the number of free param- 
eters in a model always improves the fit to the data and thus increases the value. 
However, errors introduced in the model structure at the same time may actually 
make the results of the model less reliable. For this reason, to assess the reliabil- 
ity of a model, it is necessary to ensure that the model is not over-parameterized. 

335 For this purpose, the Bayesian model probabilities provide a superior tool. To our 
knowledge, the only model that has been proven parsimonious in this sense, and 
is capable of describing litter decomposition process globally and presenting the 
results as statistical probability densities, is the one presented here. 

YassoO? is based on a biological process-oriented view of litter decomposition 

340 process although we used advanced mathematical methods to develop the model. 
A fundamental prerequisite we set to the model, was that its structure was not 
over-parametrized. This was necessary in order to determine unequivocal values 
for the parameters of the model. Such values were needed to obtain an unequivocal 
description of the biological process, test this description in the light of the litter 

345 bag data and calculate reliability estimates (probability densities) for the results of 
the model. This modelling approach we used is different from most of the earlier 
appr oaches take i i to model litter decomposit i on at extensive geographical scales 



appr oaches take i i to model litter decomposit i on at extensive geograpiiical scale 
fe.g. iBerg eraP. Il993t iMoorhead et~aD . Il999t iTrofvmow et all . 120021 : iLiski et al 



2003 : Adair et aU " 2008f K However, it was the approach that made it possible 



for us to meet the objectives of this study. Alternatively, had we chosen the so- 
called mechanistical modelling approach, it would have been impossible to calculate 
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uncertainty estimates for model predictions. The reason is that only some of the 
relevant processes can be measured directly and these measurements are only rarely 
available in global scale. 

The global patterns of foliage litter decomposi t ion ca lculated with this model, 
agreed broadly with earlier estimates ( Matthewi . 1997 ). A detailed comparison 
was impossible, because the paramet ers rep resented in the studies were somewhat 
different, soil respiration (jMatthewa . 119971 ) vs. relative effect of climate on litter 
decomposition rate (Fig. Ht- The global patterns were also qu ite similar to those of 
net primary productivity (jSchuuii . l2003HCramer et al.l . ll99 ^. This was expected 



because both processes depend strongly on temperature and precipitation. 

The results of Yasso07 were reliable enough to distinguish the differences in litter 
mass loss rates between the sites representing four climate zones (tundra, boreal 
zone, temperate zone and tropics) and two litter types (coniferous, deciduous) in- 
side each site. As a matter of fact, already a 3% difference in remainig mass values 
(relative to the initial mass) appeared as statistically significant between the tundra 
and the boreal site. This difference was caused by only a 1.6 C higher annual mean 
temperature at the boreal site. Similar climate warming is expected to take place 
as soon as over the next few decades a t the high l atitudes, and some regions have 
warmed by this much or more already (jlPCC report . 2007 ). Based on these differ- 
ences in the decomposition rates between the current litter types and the currently 
prevailing climate conditions, Yasso07 model is reliable enough for distinguishing 
small changes in the litter decomposition rates over relatively short periods of time 
in response to climate change. 

Using a large dataset of litter bag measurements and advanced mathematical 
methods of complex inverse problems we were able to develop a globally applicable 
description of foliage litter decomposition process and present this description in 
an exact mathematical form. This model, Yasso07, can be used for estimating the 
effects of changing climate and litter type on soil carbon decomposition in statistical 
sense. Our study supported the results of previous studies in that temperature 
and precipitation are the main factors affecting litter decomposition. Our results 
suggest, however, that nitrogen content of litter or site fertility parameters have a 
negligible effect on litter decomposition rates at the global scale. Because of the 
uncertainty estimates we produced for the parameter values of the model, Yasso 07 
is currently the only model that can be used for evaluating how large the changes 
in the climate and/or litter type have to be before significant statistical differences 
can be detected by modelling. 

We conclude that Yasso07 is suitable for estimating foliage litter decomposition 
of a wide variety of plant species across the global climate conditions. As a result 
of its general yet parsimonious structure, it can be used as a litter decomposition 
module in larger ecosystem models or Earth System Models. 
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Table 2. The litter bag study sites. 
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EURODECO 


32.82 


21.85 


600 


Tropical Dry 


CIDET 


55.92 


-98.62 


260 


Cool Temperate Subalpine Moist/Wet Forest 


EURODECO 


32.82 


21.85 


300 


TVopical Dry 


CIDET 


53.22 


-105.97 


476 


Cool Temperate Moist Forest 


EURODECO 


66.53 


20.18 


280 


Boreal 


CIDET 


45.92 


-77.58 


173 


Cool Temperate Moist Fbrest 


EURODECO 


60,92 


14.02 


350 


Boreal 


CIDET 


50.60 


-127.33 


100 


Cool Temperate Wet Forest 


EURODECO 


40.82 


14.48 


250 


Mediterranean 


CIDET 


54.87 


-66.65 


500 


Cool Temperate Subalpine Rain Tundra/Wet Forest 


EURODECO 


59.82 


16.55 


63 


Boreal 


CIDET 


48.03 


-123.70 




Cool T'cuiperate Wet Forest 


ELRODECO 


50.40 


13.08 


80 


Boreal 


CIDET 


51.83 


-104.92 


53G.5 


C'ool T'emperate Steppe 


ELRODECO 


58.10 


13.28 


135 


Boreal 


CIDET 


54.60 


-126.30 


1100 


Cool Temperate Subalpine Moist Forest 


EURODECO 


63.22 


14.47 


325 


Boreal 


CIDET 


60.85 


-135.20 


667 


Cool Temperate Subalpine Moist Forest 


EURODECO 


55.65 


13.32 


46 


Boreal 


EURODECO 


42.73 


-8.75 


530 


Mediterranean 
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Table 3. The properties of litter bag datasets used to calculate 

the Yasso07 model solution: litter species, numbers of measure- 
ments {Nfn) measurement timelines {At), ranges of mean annual 
temperature {Tm,min,Tm,max), amplitude of annual temperature 
variations {Ta,min,Ta,max) and precipitation {Pa,min, Pa,max), and 
initial nitrogen contents (N-cont.). Data sets with chemical com- 
position measured as a function of time are denoted by 



Name 


Species 


iV„ 


At 


r,„ ,„-.„ 


T,„ „„,-, 


r„ ,„-.„ 


r. .,„„-, 


Pa ,„.,. 


Pa ,„„,■ 


N-cont. 


Reference 








|a] 


l°C] 


l"C] 


I'C] 


I'C] 


|m] 


[m] 


|%] 




spine* 


Pinus sylvestris 


1196 


5.42 


3.8 


3-8 


11-4 


11-4 


0-722 


0.722 


0-40 


Berg ct al., 1991a, b 


Ipine* 


Pinus contorta 


128 


3.97 


3.8 


3.8 


11.4 


11.4 


0.722 


0.722 


0.40 


Berg et al., 1991a, b 


birch* 


Betula pubescens 


148 


4.00 


3.8 


3.8 


11.4 


11.4 


0.722 


0.722 


0.70 


Berg et al., 1991a, b 


liar ad* 


Pinus sylvestris 


256 


4.00 


1.3 


1.3 


12.2 


12.2 


0.650 


0.650 


0.40 


Berg et al., 1991a, b 


manja* 


Pinna sylrcilns 


240 


3-00 


1.0 


1.0 


13-1 


13-3 


0-700 


0.700 


0-40 


Berg ct al., 1991a, b 


norrl* 


Pniiia sgin ,ln3 


176 


3-33 


1-2 


1-2 


11-8 


11-8 


0-595 


0.595 


0-40 


Boig ct al-, 1991a, b 


nemie* 


Pinus sylvesLris 


236 


3.06 


6.2 


6.2 


9.8 


9.8 


0.930 


0.930 


0.40 


Berg et al., 1991a, b 


acsal 


Acer saccharum 


696 


10.06 


-7.0 


26.0 


0.8 


20.7 


0.209 


3.914 


0.81 


Gholz et al., 2000 


ambrl 


Ammophila brevUigulata 


21 


10.22 


-3.6 


22.1 


2.0 


20.7 


0.260 


3.500 


0.67 


Gholz et al., 2000 


angel 


Andropogon gerardii 


35 


9.75 


-7.0 


26.0 


2.1 


15.8 


0.284 


1.847 


0.62 


Gholz et al, 2000 


aspnl 


Populus trcmuloides 


126 


6.09 


-9.8 


9.3 


5.3 


20.5 


0.261 


1.783 


0-61 


Trofymow, 1998 


beecl 


Fagus grandifolia 


125 


6.09 


-9.8 


9.3 


5-3 


20-5 


0-261 


1.783 


0-66 


Trofymow, 1998 


belul 


Betula lutea 


24 


8.98 


-3.7 


26.0 


2-0 


17-4 


0-700 


3.500 


1-60 


Ghok ct al-, 2000 


boerl 


Bouteloua eriopoda 


35 


8.04 


-7.0 


25.6 


0-9 


20 -T 


0-233 


3.914 


0-86 


Gholz ct al-, 2000 


bogrl 


Bouteloua gracilis 


35 


8.01 


-7.0 


25.6 


0.8 


17.4 


0.209 


2.952 


0.96 


Gholz et al., 2000 


bsprl 


Picea mariana 


126 


6.09 


-9.8 


9.3 


6.3 


20.6 


0.261 


1.783 


0.70 


Ttotymow, 1998 


cedal 


Thuja plicata 


126 


6.09 


-9.8 


9.3 


5.3 


20.6 


0.261 


1.783 


0.59 


TVofymow, 1998 


cegrl 


Ccinnlhii. ,„r,,,,ii 


32 


9-00 


6.8 


26-0 


0.8 


13.8 


0.700 


2.952 


l-,33 


Gholz ct al-- 2000 


conul 


Carnii.-. inillalii 


20 


9-75 


5.5 


25-6 


0-9 


13.8 


0.310 


3.914 


0-81 


Gholz cl al-, 2000 


dflrl 


Pseudoisuga menziesii 


126 


6.09 


-9.8 


9.3 


5.3 


20.6 


0.261 


1.783 


0.65 


Trofymow, 1998 


drgn 


Drypetes glauca 


623 


10.06 


-7.0 


26.0 


0.8 


20.7 


0.209 


3.914 


1.97 


Gholz et al., 2000 


fagrl 


Fagus grandifolia 


43 


7.99 


-7.0 


25.6 


0.8 


17.4 


0.284 


2.952 


0.85 


Gholz et al., 2000 


fernl 


Pieridium aquilinum 


126 


6.09 


-9.8 


9.3 


5.3 


20.6 


0.261 


1.783 


0.82 


Trofymow, 1998 


fesci 


Fesli,,:,, hallit 


126 


6.09 


-9.8 


9.3 


5.3 


20.5 


0.261 


1.783 


0.64 


Trofymow. 1998 


jpinl 


Pinus bavksiana 


126 


6.09 


-9.8 


9.3 


5.3 


20.5 


0-261 


1.783 


1-25 


Trofymow, 1998 


komyl 


Kobresia myosuroides 


20 


10.04 


-3.6 


25.6 


0.9 


20.7 


0-233 


3.914 


1-07 


Gholz ct al-, 2000 


latrl 


Larrea tridentata 


42 


7.07 


-7.0 


26.0 


0.8 


16.5 


0-233 


3.500 


2-14 


Gholz et al-, 2000 


litul 


Liriodendron tulipifera 


25 


8.03 


-3.6 


22.1 


2.0 


20.7 


0.260 


3.500 


0.72 


Ghok et al., 2000 


Ipilb 


Pinus contorta 


623 


4.00 


0.6 


9.8 


7.7 


16.1 


0.469 


0.911 


0.40 


Berg et al., 1993 


piell 


Pinus elliottii 


187 


10.00 


-7.0 


26.0 


0.8 


20.7 


0.260 


3.914 


0.36 


Gholz et al., 2000 


pircl 


Pinna rnsinasn 


536 


10.22 


-7.0 


25-6 


0.9 


20.7 


0.209 


3.914 


0-59 


Gholz ct al-, 2000 


piall 


Pnin^ ^Irolin.^ 


25 


10.05 


-7.0 


12-7 


9.0 


16.6 


0.284 


2.291 


0-62 


Gholz ct al-, 2000 


psmel 


Pseudoisuga menzesii 


32 


7.99 


-7.0 


18.0 


9.0 


16.5 


0.284 


2.291 


0.82 


Gholz et al., 2000 


quprl 


Quercus prinus 


686 


10.22 


-7.0 


26.0 


0.8 


20.7 


0.209 


3.914 


1.03 


Gholz et al., 2000 


rhmal 


Rhododendron macrophyllum 


114 


9.03 


-7.0 


25.6 


0.8 


15.8 


0.233 


3.914 


0.42 


Gholz et al., 2000 


spall 


Spartina alterniflora 


61 


9.04 


-7.0 


26.0 


0.8 


20.7 


0.260 


3.914 


0.71 


Gholz et al., 2000 


spilb 


Pinn.s syhicslris 


672 


5.39 


-1.7 


16-7 


4-8 


15-1 


0-443 


1.500 


0-40 


Berg ct al., 1993 


tamml 


Lanx laricina 


126 


6.09 


-9.8 


9.3 


5.3 


20.6 


0.261 


1.783 


0.56 


Trofymow, 1998 


thpU 


Thuja plicata 


676 


10.22 


-7.0 


26.0 


0.8 


20.7 


0.209 


3.914 


0.62 


Gholz et al., 2000 


trael 


Triticum aestivum 


706 


10.22 


-7.0 


26.0 


0.8 


20.7 


0.209 


3.914 


0.38 


Gholz et al., 2000 


wbirl 


Betula papyrifera 


126 


6.09 


-9.8 


9.3 


6.3 


20.6 


0.261 


1.783 


0.70 


TVofymow, 1998 



